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(N ■ ABSTRACT 

o : 

O ' We report the results of a Submillimeter Array (SMA) interferometric ob- 

Tij- ■ servation of 21-/im source IRAS 07134+1005 in the CO J = 3-2 line. In order 

^ ' to determine the morpho-kinematic properties of the molecular envelope of the 

object, we constructed a model using the Shape software to model the observed 
CO map. We find that the molecular gas component of the envelopes can be 
O ' interpreted as a geometrically thick expanding torus with an expanding velocity 

■ of 8 km s~^. The inner and outer radii of the torus determined by fitting Shape 
, models are 1.2" and 3.0", respectively. The inner radius is consistent with the 

previous values determined by radiative transfer modeling of the spectral energy 
^ ■ distribution and mid-infrared imaging of the dust component. The radii and ex- 

^ ■ pansion velocity of the torus suggest that the central star has left the asymptotic 

' giant branch about 1140-1710 years ago, and that the duration of the equatorial 

■ enhanced mass loss is about 2560-3130 years. From the absence of an observed 
2 . jet, we suggest that the formation of a bipolar outfiow may lack behind in time 
^ . from the creation of the equatorial torus. 

I Subject headings: stars: AGB and post-AGB — stars: carbon — stars: imaging 

^ I — stars: individual (IRAS 07134+1005) — stars: kinematics — stars: winds, 

c3 ' outflows 



1. Introduction 

The ui iidentified 21 /xm feature was first discovered in carbon-rich proto- planetary neb- 



ulae (PPN, iKwok et al.lll989l ) and 16 21-/im sources have been found to date (IHrivnak et al. 
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20081 ). Although 20 years have past and a variety of chemical materials have been sug- 
gested as a possible carrier, no firm identification has been made. Previously propose d 



candidates include hydro genated amorp hous carbon (IBuss et al 



hydrogenated fuUerence s (IWebsterlll995l ). nanodiamonds (iHill et al.lll998l), TiC nanoclu sters 



1990 



Grishko et al.ll200l[ ). 



dSpeck fc Hofmeistei]l2004f ). etc. 



von Helden et al.l 120001) an d 0-substituted five- member carbon rings (|Papoulaiil2000l ). SiC 



Since the observation of the 21-/im feature is confined to the PPN phase of stellar evo- 
lution, one may ask whether the formation or excitation of the feature is related to the 
dynamical evolution of post-AGB evolution. For example, one may conclude that the phys- 
ical conditions (temperature and density) governing the chemical reactions may be related 
to the morphology and kinematics of the envelope. 

One of the ways to reveal the morpho-kinematic properties of PPNe is through spec- 
tral imaging of molecular lines (especially in the CO lines) of the circumstellar envelope. 
IRAS 07134+1005 (= HD 56126), a prototypical 21 /zm source, has an angular size of ~7"-8" 
(IMeixner et al.ll2004l ) and therefore is well fitted for millimeter- wave interferometric obser- 
vations. Optical Hubble Space Telescope images of IRAS 0713 4-1-1005 reveal a b right central 
star surrounded by a low surface brightness elliptical nebula (lUeta et al.ll2000l ). A near-IR 
imaging polarimetry of IRAS 07134+1005 reveals a thin (bot h geometrically a nd optically), 
limb-brightened shell wi th an equatorial density enhancenaent ( Ueta et aI1l2005 ) . In compari- 



son, t he mid-IR images (IMeixner et al.lll997l : iDayal et al.lll998l : lJura et al.ll2000l : iKwok et al. 
2OO2I ) show a hollow structure with a bright rim, with the emission region of the 21 /xm 
feature coinciding with the rim structure. 

The central sta r is variable with a period of 36.8 days, indicating a mass of 0.6 Mq 
(IBarthes et al.l I2OOOI ) . An abundance analysis by Ivan Winckel fc Reynierd (120001 ) reveals 
a metal-poor star with less than solar Fe abundance but with solar or greater than solar 
abundances of C, N, O and s-process elements, indicating the star experienced third dredge- 
up when it was on the asymptotic giant branch (AGE). The fact that it displays the family 
of aromatic infrared bands also points to its carbon-rich nature. 

In this paper, we report the results of a radio interferometric observation of IRAS 07134+1005 
in the ^^CO J = 3-2 line using the Submillimeter Arra}0 (SMA). We also present the results 
of model fitting using Shape. The outline of this paper is as follows. In Section 2 we describe 
the details of the observation and data reduction. In Section 3 we present the observational 
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results. In Section 4 we present the procedure and results of the morpho-kinematic modeling 
with Shape. In Section 5 we discuss the obtained results, and finally the present research is 
summarized in Section 6. 



2. Details of Observation and Data Reduction 

The observations of IRAS 07134+1005 were made on February 28, 2005 with SMA under 
good atmospheric conditions with a zenith optical depth at 230 GHz of 0.03-0.05. The obser- 
vations covered the frequency ranges of 344.47270-346.45638 GHz (USB) and 334.47221— 
336.45589 GHz (LSB), which include the 345.796 GHz CO J = 3-2 line in the USB. The 
array used 6 elements in the compact array configuration, with baselines ranging from 10 
to 70 m. The diameter of each antenna was 6 m, and the field of view of a single antenna 
was 34". The observation was interleaved every 20 minutes with nearby gain calibrators 
0739+016 and 0750+125 to track the phase variations over time. We calibrated the date 
using IDL-MIR, which is a data reduction software package developed by the SMA project. 
The absolute fiux calibration was determined from observations of Ganymede, and was accu- 
rate to within 15%. The final map has an accumulated on-source observing time of about 3 
hours. The single-sideband system temperature ranged from 300 to 550 K, depending on the 
atmospheric conditions and the telescope elevation. The SMA correlator had a bandwidth 
of 2 GHz with a frequency resolution of 0.812 MHz corresponding to 0.7 km s~^ in the ve- 
locity resolution at 345 GHz. The phase center of the map was taken at R.A.= 7*^16^10.26**, 
decl.= +9°59^48.0 " (J2000). Image processesing of the data was performed with the MIRIAD 



software package (jSault et al.lll995l ). The robust weighting gave a 2.28" x 1.73" CLEAN beam 
with a position angle of 2.95°. We checked the detection of the continuum emission by inte- 
grating over a 2 GHz range in LSB, but did not detect it. The upper limit of the continuum 



emission is 2.1 x 10 ^ Jy beam ^ 



3. Observational Results 



In Figure 1 we present the total intensity profile of the ^^CO J = 3-2 line. The 
profile exhibits a para bolic shape usually expected from a simple spherical outfiow (e.g., 
Morris Sz Alcocklll977l). and is about the same as those of the other CO lines (see spectra 



m, e.g. 



Zuckerman et al.lll986l : iBujarrabal et al.lll992l : iKnapp et al.lll998l . l2000l ). The peak 
intensity is 38.1 Jy at Vlsr = 74 km s~^. The velocity integrated intensity is 492.8 Jy km s~^. 
The line-width at the zero intensity level is 20 km s~^ (corresponding to an expanding veloc- 
ity of 10 km s~^). Since there is no single-dish CO J = 3-2 spectrum of IRAS 07134+1005 in 
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the literature, we cannot evaluate the extent of flux recovery by comparing with zero-spacing 
data. However, at least 70% of the flux seems to be recovered in the present interferometry. 
Since we know the size of the r nolecular envelope ( ^7"-8") measured by the BIMA obser- 
vation in the CO J = 1-0 line (IMeixner et al.ll2004l . the BIMA observation recovered 100% 
of the flux), here we assume a Gaussian brightness distribution with a FWHM of 8". Under 
this assumption 70% of the flux is recovered with our uv-coveiage of the present observation. 
Since the emission region seems to include small-scale structure, which is not resolved out 
by interferometer, the real recovery rate of the flux would be more than 70%. 

Figure 2 shows the total intensity map of the CO J = 3-2 line superimposed on the 
12.5 /xm image ta ken by the OSCIR instrument mounted on the Gemini North Telescope 
( iKwok et al.ll2002l ). The map center is taken at the position of the central star (the central 



star can be seen, for example, in the 10.3 fim and 11.7 fim images; see Figure 1 in lKwok et al. 



2OO2I ). The CO structure is clearly larger than the mid-infrared structure, and we can see a 
resolved CO structure coincide well with the 12.5 /im structure at the central region of the 
map (within ~1.5" from the map center). This ce ntral structure has a lready been hinted 
to as being hollow in the CO J = 1-0 BIMA map (iMeixner et al.ll2004l ). but the intensity 
map of Figuer 2 seems to coincide better with the mid-infrared structure than with the 
BIMA map. The central CO structure appears to be surrounded by a spherical envelope 
with an angular size of about 8". The size of the spherical envelope found in the present 
CO J = 3-2 map is almost the same as that found in the CO J = 1-0 BIMA map. If we 
assume the distance to the star is 2.4 kpc (IHony et al.ll2003l ). 8" is translated to a linear size 
of 2.9 X 10^^ cm. The spherical envelope seems to be slightly elongated in the north-south 
direction, but this elongation is due presumably to the beam pattern exhibiting elongation 
in this direction. 

Figure 3 shows velocity channel maps in the CO J = 3-2 line. The feature seen in 
Figure 3 clearly exhibits a systematic variation in velocity. The feature increases in size 
as one approaches the systemic velocity (~73 km s~^), suggesting a spherically expanding 
envelope. In addition, the central structure also exhibits a systematic variation in velocity. 
In the channels of 67-69 km s~^, an arch-like feature (convex upward) is seen in the northern 
side of the map center, in the channels of 70-73 km s~^ two intensity peaks are seen at the 
east and west sides of the map center, and in the channels of 74-78 km s~^ another arch-like 
feature (convex downward) is seen at the southern side of the map center. The p-v diagrams 
presented in Figure 4 also shows this systematic variation of the CO envelope. A reasonable 
interpretation of these variations could be an expanding torus surrounded by a spherically 
expanding envelope if the torus has a certain inclination angle. We will further consider the 
morpho-kinematic properties of the envelope in the later sections by constructing models. 
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4. Morpho-Kinematic Modeling with Shape 

In order to acquire a better understanding of tlie morplio-kinematic properties of IRAS 
07134+1005, we hav e constructed a model using the Shape software package (see e.g., 



Steffen &: L6pezll2006l ). Shape is a morpho- kinematic modeling tool used to create 3D models 
of astronomical nebulae. The modeling strategy behind Shape is to build a 3D model of the 
object, create channel maps and p-v diagrams, compare with the observational data and 
then refine the model if needed. The comparison of the generated channel maps and p-v 
diagrams with those observed is done by visual inspection. Shape does not calculate the 
numerical hydrodynamic evolution or radiation transfer equations. Since the CO J = 3-2 
line seems to be optically thick, we are only concerned with the morphology and kinematics, 
and did not try to reproduce the intensity distribution and line profile. We assumed an 
axial symmetric geometry in an effort to reduce the number of model parameters. We also 
assumed, as outlined in Section 3, that the envelope consists of an expanding torus and 
expanding sphere; both of which share the symmetry axis. 

The modeled sphere has a fixed outer-radius of 4" and a constant, radially expanding 
velocity of 9 km s^^. To model the torus, we adjusted 5 different parameters: inner and 
outer radii, thickness (height), inclination angles of the symmetry axis and the expansion 
velocity. The best-fit parameters were determined by (educated) trial and error until the 
reproduced maps closely matched the observed. The parameters obtained are: inner radius 
of 1.2", outer radius of 3.0", thickness of 4.5", and a constant radial expansion velocity of 
8 km s~^. The inclination angle was determined to be 38° to the line of sight, and rotated 
8° from the North in counterclockwise (i.e., projected position angle is 8°). In Figure 5 we 
present the polygon-mesh view of the model. The volume of the sphere and torus are each 
sampled with 20000 particles, which creates a sufficiently smooth map for comparison with 
the observations. 

In Figures 6 and 7, we present the channel maps and p-v diagrams of the best-fit model 
described above. Since Shape can currently only use a circular beam for simulating image 
convolution, we used a circular beam with a diameter of 2" instead of the actual elliptical 
beam of 2.28" x 1.73". The use of a circular beam should not affect the main features of 
the model maps, but could affect the details. The features seen in the model maps seem to 
roughly reproduce the observation, even though some noticeable differences are evident. The 
major difference is in the intensity distribution which is especially clear in the p-v diagram. 
In the observed p-v diagram (Figure 4), the intensity is clearly stronger in the red-shifted 
(high velocity) side, while the model p-v diagram exhibits a symmetric intensity distribution 
in velocity with respect to the systemic velocity (~73 km s~^). 

This asymmetry in the intensity may be due to radiative-transfer effects. The velocity 



- 6 - 



of the near side (approaching side) of the torus is, more or less, similar to that of the 
approaching side of the sphere, and therefore the emission from the near side of the torus is 
selectively absorbed by the sphere. On the contrary, the emission from the far side (receding 
side) of the torus is not absorbed by the sphere, because the velocity of the far side is clearly 
different from that of the near side of the sphere. 

Apart from the asymmetry seen in the p-v diagrams, we see further differences in the 
fine structure of the channel maps. For example, the model map shows a smooth ring-like 
feature in the 73 km s~^ channel, while the observational map clearly shows two intensity 
peaks at the east and west sides of the map center. A possible explanation for this kind 
of difference in structure could be a distorted shape to the torus. To check this possibility, 
we constructed a model with a distorted torus and found that this distortion may indeed 
explain the observed map. Of course, using a distorted model increases the number of free 
parameters and therefore has a better chance of matching the data. Nevertheless, showing 
the results of the distorted model may be valuable in paving the way for improvements to 
the simple, symmetric torus model. 

In Figure 8 we present the distorted torus model and the resulting maps. The top-right 
panel shows the geometry of the torus. A sphere with a radius of 4" similar to that of the 
symmetric model has been omitted for clarity. The torus is basically a sphere with a hollowed 
out spheroidal region. This torus has a radially expanding velocity field with a Hubble type 
law of v{r) = l/3r. The length of the major axis of the spheroid is 3.6", and the length of 
the minor axis is 1.2". The center of the spheroid is shifted to the south (along with the 
symmetry axis) by 0.7" with respect to the center of the expanding sphere. Therefore the 
one side of the hollow region is open and the other side is closed. The inclination angle of 
the symmetric axis was determined to be 55° to the line of sight, and rotated 14° from the 
north in counterclockwise (i.e., projected position angle). This distorted torus model seems 
to reasonably reproduce the fine structure seen in the observed maps. 



5. Discussion 



The firs t glimpse into t he inn er region of the molecular envelope of IRAS 07134+1005 
was made by lMeixner et al.l (120041 ) in their BIMA CO J = 1-0 observation, revealing a cavity 
structure at the central part of the envelope. In the present observation, with the higher exci- 
tation temperature and critical density of the CO J = 3-2 line, we have more clearly revealed 
the innermost structure, which is interpreted as an expanding torus. The existence of a torus 
is relatively common in evolved stars and PNe. I n fact, tori have been found in various kinds 



of objects such as AGB stars and PNe (see, e.g., iBujarrabal et al.lll998l : iHuggins et al.ll2000 
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Hirano et al.ll2004l : iHuggins et al.ll2004l : lNakashimall2005l : IChiu et al.ll2006l : iNakashima et al. 



20071 ) , which have been considered a remn ant of an equa torial enhanced mass loss developed 



in the last moment of the AGB (see, e.g., iHuggind 120071 ) . Therefore, the morpho- kinematic 



information of the torus found in IRAS 07134+1005 possibly provides us useful information 
for considering its mass-loss history in the late AGB and early post- AGB. Here, we first try 
to estimate some physical parameters of the torus of IRAS 07134+1005 by making use of 
the present data, and then we discuss the meaning of the morpho-kinematic properties of 
IRAS 07134+1005 in the scheme of the post-AGB evolution, and also discuss a possible re- 
lation between the morpho-kinematic properties of the envelope and the origin of the 21 /xm 
feature. 



5.1. Inner Radius and Lower Limit of the Torus Mass 



The inner radius of the torus is one of the most import ant pararn e ters, because it 
directly constrains the time since the central star left the AGB. IVolk et al.l (119991 ) predicted 
the inner radius of the molecular envelope by fitting radiative transfer models to the spectral 
energy distribution (SED), and gave the value of 1.2" (the original value given by Volk et 
al. is 5.9 X 10~f_pc kpc~^; here, we assume the distance of 2.4 kpc to translate the units). 



Kwok et al.l (|2002| ) gave a dust inner radius of 0.8"-1.2" from their direct imaging of the 
dust component by mid-infrared imaging. In the present research, the best-fit Shape model 
(simple, symmetric model) gives an inner radius of 1.2". This value is basicall y consistent 



with the previous values, but somewhat larger than the dust radius given by iKwok et al. 



(120021 ). A reasonable interpretation for this might be that there are no molecules in the 
region within = 1.2", even though this might be due alternately to excitation effects 
(i.e., CO molecules exist in the region i?j„ < 1.2", but are not bright in the CO J = 3-2 line 
due to high temperature and/or high density). Therefore, further CO observations in the 
higher transition lines will give us a deeper insight on this point. 

If we assume that the expanding velocity of the torus is a constant value, 8 km s~^, 
which is determined by fitting Shape model, the dynamical time scale of the inner edge 
of the torus (-Rm = 1.2") is estimated to be 1710 yr (also at the distance of 2.4 kpc). If 
we assume t he inn er radius of 0.8", which corresponds to the dust inner radius given by 
Kwok et al.l (120021 ). the dynamical time scale is estimated to be 1140 yr. On the other 
hand, the outer radius of the torus determined in our Shape modeling is 3.0", giving a 
dynamical time scale of 4270 yr under the assumption that the expanding velocity is also 
8 km s~^. (Since the outer edge of the torus seems to interact with the outer spherical 
envelope, the time scale of 4270 yr might include significant uncertainty). A straightforward 
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interpretation of these time scales is that the mass loss of the central star has gone into the 
equatorial enhanced mass-loss mode, forming a torus 4270 years ago. The phase of mass 
loss then continued for about 2560-3130 years, and finally stopped 1140-1710 years ago. 
This kinematic history of the envelope is consistent with the time-variation of a mass-loss 
rate assumed in previous radiative transfer calculations, whi ch reasonably reproduce the CO 
observations (IMeixner et al.ll2004l : iHrivnak fc Bieging]|2005l ). 

We roughly estimate the lower limit of the torus mass from the integrated intensity of 
the observed CO J = 3-2 line. (Since the CO J = 3-2 line is most likely optically thick, 
here we estimate just a lower limit of the r nass). To calcu l ate th e lower limit of the mass, 
we used the optically thin formula given by iHuggins et al.l (119961 ) . To derive the flux of the 
torus, we integrated the flux over a circle region (radius is 1.5") placed at the map center 
(the derived value is 265.6 Jy km s^^). For a distanc e of 2.4 kpc and a represen tative CO/H2 
abundance of 7.4 x 10~^ for carbon-rich envelopes (IHrivnak fc Bieging]|2005l ). we obtained 
the value Mtorus > 3.8 x 10~^ M0. If we assume that the duration of the equatorial outflow is 
2800 years as derived above, the lower limit of the mass-loss rate in this phase is calculated 
to be 1.4 X 10~^ M0 year~^. This lower limit of the mass and mass- loss rate are 2 orders 
of magnitude smalle r than values given by radiative transfer equation calculations (see, e.g., 
Meixner et al.ll2004l ). This presumably means that the CO J = 3-2 line is significantly 



optically thick. 



5.2. Meaning of the Morpho-Kinematic Properties 

A noteworthy characteristic in the morpho- kinematic properties of IRAS 07134+1005 
would be the absence of a jet (i.e., a high- velocity bipolar flow). A jet consisti ng of molecula r 



or ionized gas is often found in evolved stars and PNe together with a torus (iHugginsI 120071 ) . 
but it is not found in IRAS 07134-1-1005. The jet might be resolved out by interferometry, 
but even in such a case the existence of a high-velocity jet should be hinted in the line 
profiles given by previous single-dish observations as a high-velocity wing. However, previous 
CO observations have not found the evidence of such a high velocity wing. Therefore, 
the jet seems to be really missing in the envelope of IRAS 07134+1005, unless it is still 
not accelerated enough to form a wing component in the line profile, or fully ionized and 
slipped from the molecular line observations (i.e., too weak to be detected in molecular lines). 
Otherwise, there is no jet at any time. 



HugginsI (120071 ) recently suggested, based on statistical analyses of evolved stars and 
PNe, that a jet is formed after a torus is formed (although the formation mechanism of jets 
and tori is still not clear), and that the time lag between the torus and jet formation (here- 
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after, "jet lag" after the fashion of the Huggins's paper) is several hundred years (the median 
is 300 years but with a large scatter). The jet lags of the Huggins's sample ranged from 130 
to 1660 years. Our kinematic analyses suggests that IRAS 07134+1 005 seena s to ha ve left the 
AGB roughly 1400 years ago. According to the range of jet lags in iHuggind (120071 ). it would 
not be unnatural if a jet is well developed in the molecular envelope of IRAS 07134+1005. 
Therefore, the absence of a jet might mean that the jet lag of IRAS 07134+1005 is relatively 
long, if the jet is yet to be formed. 

Interestingly, another prototypical example of the 21-/im source, IRAS 22272+5435 
(IKwok et al.lll989l ). also shows no sign of a high- velocity jet, while this object seems to 
have a developed torus. Single-dish line profiles of the CO rotational lines exhibit an un- 
exceptional parabolic sha pe, in which no hints of a high velocity wing are found (see, e.g., 
Hrivnak fc Bieging) 120051 ). Fitting of radiative transfer models to the single dish spectra 
suggested that a high mass-loss rate (~ 10"'^ Mq) is required in the last 2000 years, sug- 
gesting that a developed torus, which is formed by the superwind, exists in the inner region 
of the molecular envelope. Although radio CO imaging did n ot resolv e the torus struc - 
ture due presumably to its sm all angular si ze fjFong et al.l 20061) . optical fjUeta et al.ll2000l ). 



near-IR (IGledhill et al.l l200lh and mid-IR flMeixner et al. 



1997; UetaetaL 



200ll ) imaging 



found a toroidal dust shell in the innermost part of the envelope. The resemblance of 
IRAS 22272+5435 to IRAS 07134+1005 tells us that all 21-/xm sources might exhibit a large 
jet lag, and also that all the 21-/im sources are at an intermediate phase between the torus 
and jet formation, even though there still could be other interpretations (i.e., the jet is/was 
very weak, or there is no jet formed at any time). 

The origin of the jet lag is not clear. However, the long jet lag potentially constrains 
the physical conditions of the innermost region of the envelope, which cannot be resolved 
by any existing equipments, if we assume a certain model. For example, one of the possible 
models, which explains the format ion of jets and tori (and jet lag), is an a ccretion disk on a 
primary /secondary star (see, e.g.. iMorrid 119871 : ISoker fc RappaportI l2000l ) . In the accretion 
disk model, the jet lag is tightly related to the viscosity of the disk: the longer the jet lag, 
the lower the visco sity [the mathematical relation between the jet lag and viscosity is given 
by iHugginsI (120071 ): see. Equation 2 in his paper]. If the viscosity is related to the chemical 
composition of the envelope, we might be able to constrain the carriers of the 21 /xm feature. 



5.3. Short Remarks on the Distorted Torus Model 



As stated in Sect 4, we have to be careful in the interpretation of the distorted torus, 
because the distorted shape of the torus increases model parameters, and also because ob- 
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servations in optically thick lines, such as CO J = 3-2, might not be the best for comparison 
with the detailed Shape model. However, the distorted torus model suggested in this paper 
seems to include an important indication. It is that the observation might be reproduced 
with the asymmetry of the torus in the axial direction, which presumably means that the 
two jet flows have not been formed simultaneously (in other words, not synchronized), while 
another interpretation still might be that the density distribution of the envelope is not ho- 
mogeneous in the axial direction (i.e., dense in the north, but thin in the south). Whatever 
the case, the torus morphology provides us useful information about the physical conditions 
of the envelope, and/or helps to constrain the formation mechanism of jets. 

6. Summciry 

This paper has reported the results of an SMA observation of IRAS 07134+1005 in the 
CO J = 3-2 line. We also performed morpho-kinematic analyses with Shape. The main 
results of this research are summarized below: 

1. The innermost structure of the molecular envelope was clearly resolved, and it is inter- 
preted as a geometrically thick expanding torus with an expansion velocity of 8 km s"^. 

2. The inner radius of the torus determined by the Shape model is 1.2", suggesting that 
the central star has left the AGB about 1400 years ago. The derived inner radius 
of the torus is consistent with the previous values determined by radiative transfer 
calculations and mid- infrared imaging of the dust component. 

3. The outer radius of the torus determined by the Shape model is 3.0", suggesting that 
the central stars has gone into the supcrwind phase from 4270 years ago. In conjunction 
with the time scale of the inner edge of the torus, the duration of the superwind is 
estimated to be 2560-3130 years. 

4. There is no jet observed in IRAS 07143+1005. The jets could be weak or absent, or 
have a long delay times, and these may be characteristic of the 21 /xm sources. 
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Fig. 2. — Total flux intensity ma p in the ^^CO J = 3—2 line superimposed on the mid- 
infrared 12.5 jjm image taken from lKwok et al.l ( 2002 ). The contour levels are 3, 19, 35, 67, 
83, 99, and 105 a, and the 1 a level corresponds to 4.80 x 10^^ Jy beam~^. The dashed 
contour correspond to —3 a. The FWHM beam size is located in the bottom right corner. 
The origin of the coordinate corresponds to the phase center. The dotted green arrows 
indicate the directions along which the p-v cuts shown in Figure 4 were taken. 
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Fig. 3. — Channel maps of the SMA compact array data in the ^^CO J = 3—2 hne. The 
velocity width of each channel is 1 km s~^ and the central velocity in km s~^ is located in the 
top left corner of each channel map. The contours start from the 3 a level, and the levels are 
spaced every 7 a. The 1 a level corresponds to 2.35 x 10~^ Jy beam~^. The dashed contour 
correspond to —3 a. The FWHM beam size is located in the bottom right corner of the last 
channel map. The origin of the coordinate corresponds to the phase center. 
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Fig. 4. — Position-velocity diagrams of tlie CO J = 3-2 line. The contour levels are 1.17, 
2.32, 3.47, 4.62, 5.77, and 6.92 Jy beam~^. The direction of the cuts is indicated in Figure 2. 
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Fig. 5. — Polygon-mesh image of the best-fit model. 2" corresponds to 7.2 x 10^^ cm at the 
distance of 2.4 kpc. 
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Fig. 6. — Channel maps of the best-fit model. 
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Fig. 7.- 



Position-velocity diagrams of the best-fit model. 
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Fig. 8. — Top left: Position-velocity diagrams of the distorted torus model. Top right: 
Polygon-mesh image of the distorted torus model (sphere is omitted). Lower panels: Channel 
maps of the distorted torus model. 



